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ABSTRACT 



CO 

■ The basic cosmological predictions of a 'polarizable-vacuum' theory of gravity 

cn 

their compliance with current observations. A Friedmann-like equation is derived 
for the theory that differs in a significant way from that of GR, for which it is 
demonstrated that a big bang is not a solution: the initial scale factor may 
be small but not zero, and must have zero first derivative. Detailed graphs of 
the predictions of the theory for combinations of the density parameters are 
given. It is shown that these predictions are compatible with current estimates 
for the age and deceleration parameter, and the maximum observable redshift, 
provided the matter density Vt m is close to the luminous matter density, i.e. 
provided there is not a significant amount of missing mass. For these best fit 
parameters, the theory predicts a maximum initial temperature that is too small 
for BBN and probably too small to support an era of radiation domination. The 
theory predicts a radiation power from a binary system that is 2/3 that predicted 
by GR, and so incompatible with observed orbital decay rate of PSR 1913 + 
16. There is some discussion of possible developments, including reconciliation 
through accommodation of vacuum anisotropy, and of the relation between PV 
and the Yilmaz theory and the quasi-steady-state cosmology. 



1. Introduction 

That the action of gravity on light can be interpreted as a variable vacuum refractive 
index has been a persistent theme in the literature - Eddington (1920); Gordon (1923); 
Balazs (1957); Skrotskii (1957); Plebanski (1959); Atkinson (1962, 1965); Bertotti (1966); 
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Winterberg (1968); Volkov, Izmest'ev & Skrotskii (1971); de Felice (1971); Evans, Nandi & 
Islam (1996a,b) - and in some of the standard texts on GR - Landau & Lifshitz (1980); Weyl 
(1960); M0ller (1962); Fock (1964). The condition that gravitational action on light be exactly 
expressible in terms of a variable scalar refractive index n(x, t) is just that the metric can be 
put into isotropic form, whereupon c/n(x, t) is the coordinate speed of light (the proper speed 
of light is always c). Since this condition is met outside the event horizon of the Schwarzschild 
solution, and by all the FLRW metrics, the refractive index interpretation and corresponding 
mathematics of EM propagation in gravitational fields has wide application. In the event 
that the coordinate system cannot be put into isotropic form, de Felice (1971); Volkov, 
Izmest'ev & Skrotskii (1971); Landau & Lifshitz (1980) have given the general relationship 
between the metric tensor and the generalized constitutive relations to be ascribed to the 
vacuum that will correctly predict the behaviour of the EM fields. These relations include the 
requirement, in order that EM radiation suffer no bi-refringence, that the electric permittivity 
and magnetic susceptibility are mutually proportional. In summary, these works are based 
upon the premise that, ignoring the action upon matter but otherwise without loss of physical 
accuracy, the gravitational action on EM fields in a given metric can be interpreted as due 
to gravitationally-modulated vacuum polarizability. 

Wilson (1921) was probably the first to suggest extension of the (specifically) refractive 
index paradigm to include matter, by which means he demonstrated agreement with theory 
to the level of the Newtonian approximation. Subsequently, Atkinson (1962, 1965) gave the 
dependence required of the rest mass on the local vacuum refractive index in order that 
matter precisely follow the geodesies of the Schwarzschild solution. Interpreted literally, 
though the details are omitted, the efforts of Wilson and Atkinson amount to the presumption 
of an electromagnetic basis for mechanical mass. Supposing this claim can be substantiated, 
such an approach, however, does not constitute an electromagnetic theory of gravity for as 
long as the Einstein action has no electromagnetic explanation. Barcelo (2001) have given 
suggestions for how the latter might emerge from an effective action wherein the observable 
degrees of freedom are really linearized fluctuations of an invisible background field. 

In addition to the authors cited above, Dicke (1957, 1961), also, was an advocate of the 
refractive index interpretation of the action of gravity on light. Like Wilson and Atkinson, 
Dicke went on to speculate on an electromagnetic basis for mechanical mass, whereby the 
polarizable vacuum paradigm could be extended to encompass not only EM but also matter. 
His work has recently been re- visited by Puthoff (2002). Unlike all previously cited authors 
apart from Puthoff, however, it is evident in those papers that Dicke's motivation was to 
demonstrate the viability of a theory that is electro magnetically motivated from the outset, 
rather than simply to forge agreement between GR and a theory of a polarizable vacuum 
(hereafter, TV). An important consequence of that difference is that Dicke adduces a metric 
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that differs from that of Schwarzschild in the higher order terms of the PPN expansion. 
Specifically, if the components of the line element for the latter are written in isotropic 
coordinates as 

ds 2 = {^^j C 2 dt 2 - (1 + rfx 2 (1) 

where = Gm/2rc 2 , then the corresponding PV line element, (Puthoff 2002), is 

ds 2 = e-^c 2 dt 2 -e^dx. 2 . (2) 

(Here and throughout, rather than set to 1 through redefinition of the time coordinate, the 
speed of light c is retained in respect of the interpretation of K as modifying the vacuum 
speed in a Euclidean background.) Consistent with this line element, the PV metric always 
has the form 

ds 2 = c 2 dt 2 /K - Kdx 2 (3) 

where K = K (x, t) is the vacuum polarizability. The specification of PV is complete by 
assigning to the If -field, on the basis of maximum mathematical simplicity consistent with 
local Lorentz Invariance, the action of a scalar wave equation with speed c/K. 

It is important that though Dicke's reasoning lead him to equations (2)) and (3), it 
did not do so unambiguously; his choice rested ultimately on the criterion of mathematical 
simplicity, given that gravitational action is fundamentally electromagnetic and mediated by 
a polarizable vacuum. This is best demonstrated by contrasting the very simple expressions 
obtained by PV for the refractive index and mass-dependency on the gravitational field 
strength parameter ifj, for example, with the results of Atkinson (1962). In the latter, it 
is recalled, GR effects on EM and matter are driven by the requirement that the (GR) 
Schwarzschild metric be reproduced exactly. Atkinson's work achieves the stated aim, but 
at the price of cumbersome expressions relating the refractive index and rest-mass to ip, 
making the correspondence look like a force- fit (Atkinson 1962). In particular, eliminating 
■0, it can be shown that in the case of the Schwarzschild metric the rest-mass is a sixth- 
order polynomial in the refractive index, and cannot, therefore, be given in closed form, 
de Felice (1971), whose aim was to establish the formal correspondence for the EM fields 
alone, concludes likewise. Specifically he writes that "the medium equivalent approach does 
not look any simpler than general relativity; on the contrary it may be mathematically 
more cumbersome and less elegant". By contrast, in PV the refractive index and mass- 
dependencies are both exponentials of ip, and therefore are related by a simple power. 

As noted by Dicke (1961), the Yilmaz theory of gravity (Yilmaz 1992; Cooperstock 
& Vollick 1996; Alley 1995) has in common with PV a line element of the form (3) in the 
particular case that T 00 is the only non-zero element of the stress energy tensor, i.e. where the 
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source matter is static, and the pressure / stress and momentum terms are zero (Alley 1995). 
Given that Yilmaz originally conceived of his theory in terms of a polarizable vacuum, Yilmaz 
(1986), this shared characteristic lends some support to the 'electromagnetic interpretation' 
of the PV metric. Further connections between the two theories are discussed in section 5. 

Though based upon a single scalar function, this remarkably simply theory passes the 
four standard tests of GR. Yet, to the extent that PV is driven fundamentally by an EM 
paradigm, it is a representative of a genre that is quite different from GR. However, apart 
from the recent report by Puthoff, it appears not to have been otherwise pursued since 
Dicke's original publication. The aim of this document is to fill the gap and thereby permit 
a more informed assessment of the viability of PV and its domain of applicability. In the 
following then, the theory is applied to a few cosmological and astrophysical problems, and 
its predictions compared with the data. 

2. PV Friedmann equation 
2.1. Action 

Dicke and Puthoff give the classical action for a single massive charge plus EM fields in 
a polarizable vacuum (in SI units) 




where x(t) is the particle trajectory, v = rfx (t) jdt the particle velocity, and v its speed. If 
in the integrations one takes dxo = cdt, then the action (equation (4) and subsequently) is c 
times its usual definition. The particle mass tuq is independent of K. The EM field strengths 
E and B have the usual relation to the potentials, E = — (V<f) + dA/dt) , B = V x A. 
The contribution from the i^-field itself is that of a classical scalar wave-equation in curved 
space-time: 

4 f 

l k = J d 4 xV=gg ab d a Kd b K (5) 

where, from equation (3), in Cartesian coordinates g ab = diag {K, —l/K, —l/K, —1/K}. Ev- 
idently this action is not invariant under arbitrary coordinate transformations, and therefore, 
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unlike GR, the background coordinates here contain some physical significance. 

Several obvious changes and generalisations of equation (4) are required to render PV 
compatible with the Cosmological Principle: 

1. Long-range EM forces are negligible compared to gravitational forces, whereupon the 
EM interaction part of the action due to any charge can be safely ignored. 

2. The single charged particle must be replaced by a uniform distribution of neutral mass 
that is static in the background frame, whence v — 0. 

3. moc 2 S 3 (r — r (t)) — > p m where p m is a uniform (in the background frame) and constant 
energy density of matter. Subsequently we will choose K — 1 at the present time, so 
that p m will be the present energy density of matter. 

4. On the cosmological scale, the field K is assumed to be isotropic, homogeneous, and 
therefore a function of time only, for which it will be convenient to make the substitution 
K = a 2 (t), which amounts to the selection of the flat-space metric out of the three 
FLRW space-times. The background to this choice is discussed further below. 

5. In anticipation of the need to accommodate accelerated expansion, as in conventional 
GR the action will be supplemented with a cosmological (henceforth 'vacuum') term 



where p v has the units of energy per unit volume (the energy density of the vacuum 
when a = 1). 

With these changes the PV action, (4), becomes 




(6) 




- a 2 e E 2 



(7) 



and consistent with equation (3) the line element is 

ds 2 = c 2 dt 2 /a 2 (t) - a 2 (t) dx 2 . 




2.2. 



Hamiltonian 



Carrying out the space integrations, the Lagrangian is 

t \r\ a2 ° 2 ( da \ 2 2 Pm , T \ 



(9) 
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where the electromagnetic part is 

Lem = ^j f d 3 x L 2 e E 2 -4—)- ( 10 ) 

The corresponding contributions to the Hamiltonian are 

.dL a 2 c 2 ( da\ 2 2 p m TT 

H = a—-L = — ( — )+p v a 2 + ^ + H EM (a) 11 
da 8nG \ dt J a 

where Hem (a) is the electromagnetic energy density. Hem (a) has the usual (i.e. GR) 
dependence on the scale factor p r /a 2 say, as may be deduced by standard arguments (Misner, 
Thorne & Wheeler 2000). With this in equation (11), the total energy density, p s say, is 

a 2 c 2 (da\ 2 2 p m p r 

' — ) +p v a + — + — = p E . (12) 



87rG V dt } ' a a 2 



2.3. Relation to the GR Priedmann equation 

In accord with equation (8) the proper time r is given by dr = dt/a (t), whereupon one 
recovers the FLRW line element 

ds 2 = c 2 dr 2 - a 2 (r) rfx 2 . (13) 

With this, equation (12) written in proper time is the familiar Friedmann equation of GR 
except for a missing factor of 3 and some changes in sign. The two equations can be united 
in the form 

3c 2 / da\ 2 2 p rn p r 

x ^g\Tt) =pva + (14) 

where in GR A = 1, and in PV A = —1/3. ps has been defined so that when p v = in GR, 
a positive value gives rise to collapse, I.E. it has the sign of in the traditional notation. 
In PV a positive p s connotes a positive net overall energy. In both cases - GR and PV - 
the magnitude and sign of both ps and p v can be decided, for example, by regression of the 
data for a versus r onto the corresponding Friedmann equation (i.e. equation (14) for either 
A). Both the mass energy density and the electromagnetic energy density could, in principle, 
be determined locally, without reference to either the Friedmann equation or the theoretical 
framework in which it is to be employed (GR or PV). In any case, the sign of both will be 
positive. 



Carrying out the usual normalisation using the present value of the Hubble constant 

1 da 
a dr 



da(r ) (15) 
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(the scale factor at the present time is 1), equation (12) can be written 

vM' = * a ' + ^ + %- nE - (l6) 

where 



8nGp 

and 



= 3c2H2 ; i e {m, r, v, S} (17) 



n v + n m + a r - fi E = a (is) 

are true at all times. Of course equation (18) can be used to eliminate any one of the 
4 fi's in equation (16), restoring it to an equation that is explicitly dependent on just 4 
observationally determined constants (3 f2's and H ). With the additional definition 

tt = Q v + VL m + VL r (19) 

employed in the (GR) literature, one has (exclusively in GR) fi s = — 1, whereupon 
Vl = 1 corresponds to flat space. But VL = 1 with f2 defined in equation (19) has no such 
special significance in PV, and therefore the substitution (19) is less useful when PV and 
GR are being compared in parallel, and so will not be employed here. Also, in the flat-space 
version of PV presently under consideration, f2s 7^ signifies non-zero overall energy, with 
no contribution from the spatial variation of K. By contrast, in the Friedmann equation 
of GR, fis 7^ always signifies a contribution from the (non-flat) geometry - it is not an 
arbitrary constant of integration. 



3. Cosmology 
3.1. Initial condition 

From equation (16) the current age as a function of the Q is 




(20) 



a 2 Vt v + Vt m /a + VL r /a? — £7s 



In the case of GR (A = 1) the operand of the square root is always positive as a — > 
because fl r is positive. Therefore in GR one is free to suppose that a (0) = was the initial 
state of affairs, and take the age as the proper time elapsed since then. By contrast in PV 
(A = —1/3) the scale factor cannot shrink to zero else the operand of the square root will 
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go negative. Denning the age as the proper time since the most recent minimum (including 
the case that the scale factor oscillates), a (0) is the smallest real positive root of 

a 2 (0) Q v + tt m /a (0) + Vt r /a 2 (0) - fi s = . (21) 

Evidently da/dr = at that point, and therefore the expansion is initially exclusively 
quadratic - i.e. is not a 'bang'. 

3.2. Very early history 

The early history is decided by the magnitude of the initial (and therefore minimum) 
scale factor given by the solution of equation (21), for which some estimate of the value of the 
coefficients is necessary. The CMB mass-equivalent energy density is 2 x 1(T 31 kg/m 3 (the 
contribution from stellar radiation is around 20% of this value and can be ignored relative 
to the imprecision of this analysis). Taking Hq ps 70 km/s/Mpc, giving a 'Hubble time' 
1/Hq ps 14 Gyr, gives (from equation (17)) Q r ~ 2 x 10~ 5 . The widest possible range for 
the matter density is 0.005 < Q m < 0.3 (absorbing therein any uncertainty in H ). The 
lower limit corresponds to luminous baryonic matter (only). The upper limit is supplied 
by, for instance, fitting super-nova data to the (GR) Friedmann equation - see for example 
Perlmutter et al (1999) and Riess (2000). Although the applicability to PV of such a large 
value for the upper limit is questionable, it will turn out not to make any difference to the 
conclusion of this section. 

To achieve a lepton era requires a temperature of around Ti epton ps 10 10 K and therefore 
an initial scale factor of no more than a (0) = TcMB/Ti ept0 n ~3x 10~ 10 . At those temper- 
atures, the matter term in equation (21) is certainly negligible compared to the radiation 
term, and therefore one can take 

a 2 (0) a + a r /a 2 (0) - ft s = . (22) 

Equation (18) gives that the magnitudes of Q v and can differ at most by a number of order 
unity, from which it is deduced that the first term in equation (22) is negligible compared 
with the second term. Therefore f2s ~ Q r /a(0) 2 ~ 2 x 10 14 . But such a large number has 
un-physical consequences on the age, for example, which can be shown as follows. With 
reference to equation (20), over the interval a (0) ^ a ^ 1, the very large Q r makes the 
square root very small everywhere except when a is very close to a(0). Since Q v and 
have similar magnitudes, it follows that a 2 \ fl v \ <^ and aVL m Q r over the region of 
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importance. Therefore equation (20) is well approximated by 

1 



1 



da 



1 



V3H J y/ttx - tt r /a 2 



(23) 



o(0) 

Rearranging a little and then performing the integral: 

i 



da 



a 



o(0) 



a/a(0)f 



« 2 (o)/(«/«(o)) ; 



(24) 



o(0) 



a(0) 
\/3f2 r i/o 



Using the numbers above, this gives an age of only 540 yrs. Consequently it must be 
concluded that PV cosmology is incapable of accommodating either a lepton or a baryon 
era. 



3.3. Radiation era 

Existence of a radiation era requires that the initial radiation energy density is at least 
equal to that of matter: a (0) = Q r /Q m . Using this in equation (21) (where Q v comes from 
equation (18), but in any case a 2 (0) \Q V \ <C and so can be neglected) gives 

fi s « 2Vt r /a 2 (0) = 2ttl l /tt r < 2.5 , (25) 

the upper limit corresponding to Q m = 0.005. The approximation used for the lepton era, 
(24), gives 

~ v^^o V 3 n m H v ; 

which, using the numbers above, gives an age less than 7.2 Gyr. However, that approximation 
required f2s ^> 1, so this result cannot be expected to be very accurate. In fact an accurate 
numerical evaluation of equation (20) with fi^ — 2.5 gives an age of 6.6 Gyr or less, the 
upper limit corresponding to Vt m = 0.005. This age is unacceptably short - too small to 
accommodate main-sequence behaviour in globular clusters for example (Peacock 1999). It 
must be concluded that PV cosmology cannot encompass an era of radiation domination. In 
short, combining the conclusions of this and the previous section, PV cosmology starts with 
a big whimper. 



3.4. Mass density 



Current estimates of the matter density based upon observation of the proportions of 
the light elements inferred from big-bang nucleosynthesis (BBN) are around Q m ~ 0.1. But 
given the absence of a baryonic era in PV, the theory is denied the chain of inference that 
leads from BBN to this value. Further, though recent observations of super-nova lead to 
£l v « 0.7, VL m w 0.3, =>- VLy, ~ - the difference between this value of Vl m and the value 
inferred from BBN is missing matter - these super-nova are due (effectively) to regression 
of the data onto the curve predicted by equation (16) when A = 1, which results cannot be 
carried over into PV, where A = —1/3. With the bases for large contributions from dark and 
missing matter gone, it seems that at least for the present PV should be examined with an 
Q m much closer to the observed - luminous - density, plus whatever component is found to 
come from the neutrino. This conclusion is vulnerable to revision however, if a significant 
dark-matter contribution turns out to be demanded for other reasons, such as, for example, 
anomalous rotation of spiral galaxies. 



3.5. Cosmological red shift 

A consequence of the non-zero initial scale factor is the upper limit it imposes on the 
cosmological redshift of distant objects: max (z) = l/a(0) — 1, where a (0) is the smallest 
positive root of equation (21) (a (0) - and therefore max (z) - are functions of the f2;). Fig. 1 
is a plot of the maximum redshift versus matter density for f2 s = and f2 s = ±0.1 obtained 
from a numerical solution of equation (21). It is observed that PV can accommodate all the 
observed redshifts provided there is an overall deficit in the energy density. At the time of 
writing the largest claimed redshift for a radio galaxy is TN J0924-2201 at z = 5.19 (Breugel 
2000). But Fig. 2 shows that if, axiomatically, (say) f2s = 0, then PV is incompatible 
with a cosmological origin of redshifts beyond 3, since these correspond to Q m < 0.005 - 
the luminous mass density. Since the observational upper limit on the redshift is subject 
to constant review, it is not possible to use that to deduce a firm relation between £7^ and 
Q m in PV. Loose constraints on the range of these parameters can be inferred from Figs. 2 
and 3, which shows simultaneously how they affect the maximum redshift, the current age, 
and the deceleration parameter over different ranges. In both figures redshifts beyond 10 are 
increasingly crowded and have been suppressed. Broadly, one sees PV requires that 

0.1 < tt s < 0.4, 0.005 < tt m < 0.1 (27) 

the intersection of these two domains giving a reasonable range of ages, and a range of 
redshifts that encompass the most distant observations to date. Note however that in light 
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of the arguments of the previous section (absence of motivation for dark or missing matter), 
current observations constrain Q m to the lower end of this range. 



3.6. Deceleration parameter 

It would appear from the form of equation (16) that the deceleration parameter q = 
—ad/ a 2 , since it is independent of the scale of both a and r, is independent of A. However, 
taking into account that the four Qi on the right hand side of equation (16) are internally 
related in a way that depends on A through equation (18), the PV prediction for the decel- 
eration parameter as a function of the three remaining true degrees of freedom differs from 
that of GR. Differentiating both sides of equation (16) one obtains 

2Xdad 2 a da ( Q m 2Qem^ 



2a n v - ^ - (28) 



and therefore 



H 2 dr dr 2 dr\ v a 2 a 3 



ad a (aQ v — Q m /2a 2 — Q r /a 3 ) 

Q -9 



a 2 a 2 Vt v + Vt m /a + VL r /a 2 — fi s 

a 4 Q v — a 2 Qj] + aQ m + Q r 
It may be observed that in PV, because the first derivative of the scale factor is zero when 
t = (it is infinite in the FLRW cosmologies), the deceleration parameter defined above is 
negative infinite there. The present (and therefore 'local') value for both theories is 

. i^^yj J 2 I i~^^> 

Qo = q\ a =i = 75 q 775 T q 

i'u i'S t i'm T i'r (30) 

= \ (Q m /2 + tt r - Q v ) = \ (3fi m /2 + 2tt r - Q s ) - 1 
A A 

where equation (18) was used. For PV in particular 

q = -3(tt m /2 + tt r -tt v ) 

= -3 (3ft m /2 + 2Q r - fi E ) - 1 . 

Rearranging equation (31), one sees that the condition for a locally accelerating universe 
(q < ) is that 

fi s < 3tt m /2 + 2Vt r + 1/3 . (32) 

From equation (27) it is observed that the plausible range of f2s and fl m encompasses both 
local acceleration and deceleration. Specifically, the range of q consistent with equation (27) 
is 

-1.15 ^ g < 0.18. (33) 
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Contours of constant go versus f2s and Vt m are plotted in Figs. 2 and 3. 

Some caution is necessary in comparing published results of q with the value predicted 
by equation (31): where q is deduced effectively from a parametric fit of the Friedmann 
equation to the scale factor as deduced from redshift and brightness data, e.g. Riess (2000), 
the reported value will be applicable exclusively to GR and not PV. 



3.7. Future collapse and oscillation 

The condition for future collapse is that the right hand side of equation (16) becomes 
zero at some a > 1, i.e. equation (21) possesses more than one real positive solution. At 
that point (da/dr) 2 = and the expansion stalls and will subsequently reverse. Combining 
equations (18) and (21) gives that future collapse requires 



a 2 tt v - 1/3 - ft v - tt m - tt r + — + -- = (34) 

a a 2 



to have a solution in the range a > 1. Solving equation (34) for Q v , 

3 (cr — 1) a (a + 1) cr 
one sees that the range 1 < a < oo corresponds to 

< n v < oo . (36) 

An infinite vacuum density corresponds to a turn-around point at the present time when 
a = 1. (A turn- around earlier than this has effectively been excluded by the normalisation 
process in which the Q defined in equation (17) are based upon on a = 1 present-day values.) 
The condition Q v = is a special 'critical' case and is discussed separately below. Inserting 
equation (36) into equation (18), the constraint on the total energy for future collapse is 

l/3 + Q m + Q r < Q s < oo. (37) 

Solutions of equation (16) satisfying equation (37) are necessarily oscillating and symmetric 
about the times of maximum and minimum scale factor. If the constraint equation (37) 
is not satisfied, then the expansion is monotonic ('open'). The locus of this constraint is 
shown as a heavy line in Fig. 3. (The line may be regarded as a contour of constant period 
for the particular case that the period is infinite - other equal-period contours could be 
drawn above this line.) It may be observed that a limited range of plausible combinations 
of age and maximum redshift lie above this line - the majority of plausible combinations lie 
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Fig. 2. — For values of Q m close to the density of luminous matter, contours showing the 
dependency on fls and Q m of: the current age in Gyr (solid ellipse-like contours), the maxi- 
mum possible cosmological redshift (broken lines), and the current value of the deceleration 
parameter q (dotted lines). Redshifts beyond 10 are not shown. 
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Fig. 3. — For a large range of f2 m , contours showing the dependency on f2 s and fl m of: the 
current age in Gyr (solid ellipse-like contours), the maximum possible cosmological redshift 
(broken lines), the current value of the deceleration parameter q (dotted lines), and the 
boundary between oscillating and open cosmologies (heavy line). Redshifts beyond 10 and 
ages beyond 25 Gyr are not shown. 
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below. Comparison of equation (37) with equation (32) reveals that monotonic expansion 
is a sufficient but not necessary condition for local acceleration. (Consequently, a local 
acceleration does not guarantee continued expansion.) 

As an example, PV predicts an oscillating cosmology using values Q m = 0.01 and 
VLy, = 0.4, and therefore - from equation (18) - Q v — 0.06. The first positive root of equation 
(21) is a (0) = 0.027, implying a maximum observable redshift of 36. Using this value for 
the lower limit of the integration in equation (20), and using 1/H — 14 Gyr, one obtains a 
present age of 14 Gyr. From equation (31), the present deceleration parameter is 0.16. From 
further numerical investigation it is deduced that for these values the expansion started to 
decelerate at an age of 7 Gyr, i.e. 7 Gyr ago. Figs. 4 and 5 show the evolution of the scale 
factor over time from numerical solution of equation (16) for these particular values. The 
oscillation has a period of 109 Gyr and a maximum scale-factor of 2.65. Also included for 
comparison is a plot of the scale factor according to GR for the currently favoured values of 
VL m = 0.3 and fls = 0, whereupon equation (18) gives that Q v = 0.7. 

3.8. Monotonic expansion 

Any combination of f2 s and Q m below the heavy line in Fig. 3 gives rise to monotonic 
expansion. If the mass density is left unchanged from the previous example at Q m = 0.01, 
and f2s is reduced by 0.1 from Q s = 0.4 to fi s = 0.3, then, from equation (18), Q v = 
—0.04. Recalling that the condition for oscillation is that the vacuum density is positive, 
this combination must lead to an open universe. The first positive root of equation (21) is 
a (0) = 0.035, implying a maximum observable redshift of 28. Using this value for the lower 
limit of the integration in equation (20), and using 1/H — 14 Gyr, one obtains a present age 
of 15 Gyr. From equation (31), the present deceleration parameter is -0.15 - i.e. a positive 
acceleration. 



3.9. Critical expansion 

If there is no contribution from the vacuum at all, then from equation (35) Q v = =>- 
a = oo , i.e. PV predicts an oscillating universe with an infinite maximum scale factor and 
therefore with an infinite period. Constructing an example using the same mass-density as 
in the previous example, Vt m = 0.01 , equation (18) gives that Vt^ = 0.34. The first positive 
root of equation (21) is a (0) = 0.031 , implying a maximum observable redshift of 31. Using 
this value for the lower limit of the integration in (19), and using 1/H — 14 Gyr, one obtains 
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Fig. 4. — Evolution of scale factor over 150 Gyr. The solid lines are predictions of PV 
(A = —1/3) given 1/H$ = 14 Gyr and assuming a mass density close to the visible value, 
Q m = 0.01, for three different values of vacuum energy density, and consequent value for 
Q-£. The dotted line is the prediction of GR (A = 1) for the currently favoured values 
Q m = 0.3, Q v = 0.7, fi s = 0. 
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Fig. 5. — Evolution of scale factor to 15 Gyr for the same set of parameters as in Fig. 4. 
The solid lines are predictions of PV (A = —1/3) given 1/Hq — 14 Gyr and assuming a mass 
density close to the visible value, Q m = 0.01, for three different values of vacuum energy 
density, and consequent value for The dotted line is the prediction of GR (A = 1) for 
the currently favoured values Vt m = 0.3, Q v = 0.7, fis = . 
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a present age of close to 14.5 Gyr. From equation (31), the present deceleration parameter is 
-0.015 - i.e. a small positive acceleration. The scale factor versus time for this configuration 
is given in Figs. 4 and 5 along with the previous open and oscillating universe configurations, 
and the currently favoured GR evolution. 

3.10. Comparisons 

For the particular values considered above, there is little difference at the present time 
between the oscillating, monotonic, and critical expansions, and the currently favoured GR 
prediction. The present-day values and slopes are within a few percent of each other, and 
the second derivative of the scale factor (and therefore qo) is close to zero for all four tra- 
jectories. It is not significant that, of the three PV scenarios considered, the oscillating 
evolution is the closest to the GR prediction, because there is sufficient leeway in the mass 
and vacuum density parameters to change the ranking. But it is noteworthy that with little 
or no vacuum term, and no clear need for 'missing' mass, PV comes close to the favoured 
GR-determined trajectory, achieving accord with current estimates of the Hubble constant, 
model-independent estimate of age, luminous mass-density, and accommodation of the max- 
imum observed redshift. 



3.11. Summary of cosmological findings 

The principal characteristics of the scalar, flat-space, PV cosmology investigated above 
are summarized below: 

1. Compared with GR, the scale factor evolves according to a Friedmann-like equation 
in which the d 2 term is multiplied by —1/3, with all the other terms appearing in the 
same way. 

2. A consequence is that PV cosmology has no big bang, but starts smoothly from a 
non-zero scale factor. 

3. The initial scale factor is too large for a baryon, lepton, or radiation era. 

4. Consequently PV can explain neither present proportions of light elements nor the 
cosmic microwave background as the result of high temperatures in the very early 
history. 

5. The non-zero scale factor gives rise to a maximum possible redshift. 
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6. PV can predict an age, deceleration parameter, and maximum possible redshift for par- 
ticular values of mass density and overall system energy compatible with observation, 
without missing mass, and without a vacuum term. 

7. PV cosmologies compatible with observation can be monotonic, critical, or oscillating. 

8. Compatibility demands that PV contain a significant contribution from the vacuum, 
or an offset to the total energy, or both. 



4. Gravitational radiation 

4.1. Gravitational potential in the far-field 

-fT-field radiation induced by the motion of 'local' matter can be decoupled from Cosmo- 
logically induced variations provided the (retarded) distance to the radiation source is small 
on the Cosmological scale. Then, neglecting any EM contribution (including pressure), and 
the vacuum term, the appropriate action is 

L - c i **{&g ( c2 < vi ° sA '» 2 - (f ) 2 ) < 38 > 

where p m (x, t) is presumed known distribution of matter and v (x, t) its local velocity field. 
Writing K (x, t) = 1 + 20 (x, t) (<f) is the dimensionless Newtonian gravitational potential) 
and going immediately to the weak- field limit \<p\ <C 1, the action is 

L = -c 




(39) 

+ o(0 2 ) x ' 



- v 2 /c 2 

To first order, the Euler equation for is a scalar wave equation sourced by the matter: 

In the case of relativistic speeds, it is easily seen from equation (39) that the right hand 
side of the above has an additional coefficient (1 + v 2 /c 2 ) j y/l — v 2 /c 2 . From application of 
the virial theorem Dicke (1961) interprets this term as the total energy. His interpretation 
rests, however, upon the assumption of particles interacting principally electromagnetically, 
whereas the 'particle' sources of mass contemplated here are stars interacting gravitationally. 
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Further, in the latter case relativistic orbital speeds are possible only in the presence of strong 
fields - the weak-field limit is necessarily non-relativistic. Consequently equation (40) cannot 
be extended to the relativistic domain, which obviates interpretation of the relativistic terms 
and discussion of the applicability of the virial theorem to this case. 

The attention here is exclusively on the power lost to gravitational radiation, for which it 
is sufficient to take the far-field limit of the retarded Green function, whereupon the solution 
of equation (40) is 



0(r,t) 
for large r. Expanding, 



|r - r'l 



(r, t) = ^- f dh'p m (r', t - |r - r'| /c) (41) 
err J 



, rV r' 2 - (f.r') 2 "\ (42) 
r 2r 2 



the far-field potential is 



^(,,,) = |/„v{l + »yl| + + o(i)}p m (r<, t-r/c) . (43) 

Since in the Laurent expansion of the gravitational potential only terms proportional to 1/r 
contribute to the radiated power, the higher order terms can be ignored. The first two terms 
of equation (43) are the monopole and dipole contributions to the potential, proportional 
respectively to the total energy and momentum of the source matter. For a closed system 
(ignoring the radiation itself), these are constant in time and cannot contribute to radiation 
since the latter is driven by non- vanishing time-derivatives of the potential (see below). For 
radiation, therefore, it is sufficient to take 



'had 



(M) = I d 3 r' (f.r') 2 Pm (r', t - r/c) . (44) 



c 4 r 



4.2. Radiation and radiated power 

The Euler equation for distant test particles gives a force proportional to V0 which, in 
the far field, is directed away from the source - the radiation is longitudinally polarized. For 
maximum sensitivity therefore, PV radiation detectors should be aligned correspondingly. 
Since this is orthogonal to the direction predicted by GR, direct detection of gravitational 
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radiation from a Cosmological source would be a decisive test in favor of one of the two 
theories. 

Disregarding the polarization, the total power radiated is the rate of energy crossing a 
closed surface far from the source. Corresponding to the Lagrangian density C, L = J d A xC , 
the energy density in the radiation field is 

dC d<f) rad 



n 



■ > rad) 



d (d(p rad /dt) dt 



c 4 
8ttG 



(v</wr + 



2,1/ u( Prad 



c 2 V dt 



(45) 



At sufficiently great distances from the source, the field momentum is exclusively radial and 
therefore the power P, say, is 



p = cdf)dn r r 2 n (<p rad ) 



c5 wW ad ) 2 + V^ x2 



(46) 



With reference to equation (44), the spatial gradient of <fi will introduce additional powers 
of 1/r when it operates on the f, and therefore does not contribute to the power through a 
surface at a large distance from the source. The spatial gradient operating on the retarded 
time does contribute however: 

(vf(t - r/c)f ^(m_LM_y. (47 ) 

Using this, and the expression for the potential equation (44), the radiated power, (46), is 

F =4^#^( fT ' If » 2 (48) 
where I is the matrix of retarded quadrupole moments: 

Iij = J d 3 r'r' l r' :j p m (r / ,t-r/c) . (49) 

The orientation average in equation (48) can be computed as follows. Since I and i are 
symmetric they can be diagonalised by a similarity transform, i = U T AU, where A is 
diagonal and U is a real orthogonal matrix representing a proper rotation. It follows that 
the integration in equation (48) can be performed in the rotated coordinate system, r — > Ur, 
with j replaced by A: 



G 



^<(r T Arr), (50) 
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where () denotes the orientation average. (The Jacobian for the transformation is det (U r ) = 
det (U) = 1 .) Writing R = f f T = R T , one has 

(f T Af) 2 =tr((RA) 2 ) 

= £ RijJ^jkRkl^li = £ Rij 2 ^jj-^H ■ 
ijkl ij 

In a Cartesian basis (Rij 2 ) has only two distinct values: one for all the diagonal terms, and 
one for all the off-diagonal terms. Therefore 

,2\ 



(52) 



= (Raa 2 - R ah 2 ) tr (A 2 ) + (R ab 2 ) tr 2 (A) 

where a ^ b. One finds by explicit calculation that (R aa 2 ) = 1/5 and (R a 2 ) = 1/15 . Using 
this and that the trace is insensitive to rotations, one finally obtains 

P= loi( 2tr C I ' 2 ) +tr2( ' I ' ) ) • (53) 
By contrast, the GR result is (Landau & Lifshitz 1980; Misner, Thorne & Wheeler 2000): 

Pgr = £tr(V) (54) 



5c 5 

where ^ is the reduced quadrupole moment tensor = I — 1 tr (I) /3 . Using that 

2 



:2 



tr ( f ) = tr | ( I — — tr ( I 



1 



2 _ 2y m 1 2( jA , (55) 



= tr( i^--itr(l) + -ltr 2 (l))' 



tr(l 2 )-itr 2 (l) 



the GR and PV results can be combined into a single expression 

G 



p (A) = eb^ ( 3 (3 + A) tr ( r ) + 2 (1 " 3A) tr ' (56) 

where, as before, the values A = 1,-1/3 correspond to GR and PV respectively. 
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4.3. Orbital decay of a binary system 

The components of the quadrupole moment tensor for a pair of objects in mutual ellip- 
tical orbit in the x, y plane are (Peters & Mathews 1963) 

1 = -/3(l + ecos^) 2 M (57) 

where M is a 3x3 matrix with components 



M xx = — 3e sin tp cos 2 ip — 2 sin 2ip 

M yy = esinip (l + 3cos 2 ^) + 2sin2^ 

M xy = M yx = 2 cos 2ip + e cos ip (3 cos 2 ip — l) 

M,; z = M zi = 



(58) 



where i G {x, y, 2} , and 

2 _ 4G 3 m 2 m 2 (mi + m 2 ) 
er (1 — e l ) 

and where the are the respective masses, ip is the phase, a is the semi-major axis, and e 
is the eccentricity of the orbits. Use of equations (57) and (58) gives 

tr(l' 2 ) =/5 2 (l + ecos^) 4 (e 2 sin 2 ^ + 8(l + ecos^) 2 ) (60) 

and 

tr 2 (l) =/5 2 (l + ecos^) 4 e 2 sin 2 ^. (61) 
Substitution of these into equation (56) gives that the total power is 

P (A, VO = ^ (1 + e cos ^) 4 X (A, ^) (62) 

where 

X = (ll-3A)e 2 sin 2 ^ + 24(3 + A)(l+ecos^) 2 . (63) 
Since interest is restricted to secular effects, it is safe to work with the power averaged over 

r 

an orbital cycle, r, say: P (A) = f dtP(\,ip) jr. Using that, for any function /, 



T 2tt 

j dtf (*!>) = j <WW>)/^ (64) 



and that the rate of change of phase of an elliptical orbit is 

ip = k (1 + ecosip) 2 (65) 
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where k is a constant, the time-average of equation (62) is 



2tt 



J#(l + ecosV0 ^(^,-0) 

? w = pS • (66) 



/ (1 + ecos?/>) 



-2 







Performing the integrations and recalling equation (59) one obtains 

P(A) = f g^r%i±^) /(A,e) (67) 

where 

(227 + 69A)e 4 + 4(443 + 141A)e 2 + 192(3 + A) 
768 (1 - e 2 ) 7/2 

is the model-dependent 'enhancement factor', with coefficient chosen so that /(1,0) = 1. 
It is easily verified that the familiar (GR) result of Peters and Mathews is recovered when 
A = 1. The ratio of the PV to GR radiated powers, 

P (-1/3) / (-1/3, e) 1 /51e 4 + 396e 2 + 128 , 

(69) 



P(l) f(l,e) 2 V 37e 4 + 292e 2 + 96 

is plotted in Fig. 6 as a function of the eccentricity. The value of the ratio in the limit of 
perfectly circular orbits (e = 0) is exactly 2/3, and changes only by 1.5% across the full 
range of eccentricity. That the ratio is less than one is to be expected given the reduced 
number of polarizations available to the vacuum if-field compared with the two independent 
polarizations of weak-field GR radiation. This is in accord with divergence between GR and 
other theories that have different rules governing generation of radiation (Will 1977). 



4.4. Test against observations of orbital decay of a binary system 

Using the Kepler results that the energy of the binary system is E — —Gm 1 m 2 /2a, and 
its period is r b = 2ira 3 / 2 j \/G (m 1 + m 2 ), one has r b /r b = —3P/2E and therefore (Damour 
& Taylor 1991; Wagoner 1975) 

192tt (2nG(m 1 + m 2 )\ b/ \ . 
Tb = -—»{ V feC 3 ) /(A ' 6) (70) 

where v = mim 2 / (m 1 + m 2 ) 2 , and where the period has been used to eliminate the semi- 
major axis. Estimates for the masses may be obtained from the rate of periastron advance, 
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u, and the time-dilation parameter 7 in the case there is apsidal motion (Blandford & 
Teukolsky 1976). It is important that though inference of these quantities from observation 
is potentially model dependent, they are computed using the static metric of the center of 
mass system, in which case the PPN expansion of PV is known to give the same result 
as GR to the required order (Puthoff 2002). Consequently, the PV-dependent observational 
inference of these quantities, and therefore mi and m 2 , is the same as that of GR. Then, if one 
regards (66) as a theoretical prediction for the decay rate f b in terms of the observationally 
determined quantities r b , mi and m 2 , it is concluded that the prediction from PV differs 
from that by GR by the factor given in equation (69) of approximately 2/3. But since the 
observed decay rate of PSR 1913 + 16 agrees with the prediction of GR to better than 
1%, with more recent research concerned with checking relativistic corrections and other 
refinements of the analysis presented here (Taylor & Weisberg 1989), it follows that PV is 
incompatible with observation. 

4.5. Summary of gravitational radiation findings 

PV predicts gravitational radiation with the following properties 

1. Radiation is polarized along the direction of propagation; there is just one polarization 
degree of freedom. 

2. The orbital decay due to radiation of a binary system is always close to 2/3 that 
predicted by GR, with the exact relation as a function of eccentricity given by equation 
(69), and is plotted in Fig. 6. 

3. This prediction is at variance with the observational data for PSR 1913 + 16. 

5. Developments 

5.1. Anisotropic polarizability and the Yilmaz theory 

At the time of Dicke's writing, his theory passed the then four standard tests of GR. 
Dicke did not investigate the cosmological predictions of the theory, which are reported here 
for the first time. The major shortcoming of PV reported immediately above - its failure 
to predict the correct (observed) orbital decay rate of a binary system due to radiation 
- was unknown to Dicke. The factor of 2/3 disagreement between PV and GR possibly 
can be attributed to the reduced number of polarization states permitted to the radiation. 
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Fig. 6. — Ratio, PV/GR, of gravitational radiation enhancement factors. The original 
definition by Peters and Mathews is the ratio of radiation power from a binary system with 
elliptical orbits to that from a system with circular orbits. 
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Depending on the details, a theory accommodating vacuum anisotropy has the potential 
to accommodate two independent transverse polarization degrees of freedom of radiation, 
thereby bringing about agreement with the binary pulsar data (and GR). In fact Dicke briefly 
considered the possibility of an anisotropic version of the theory as plausible, but argued in 
favor of the scalar version, since it was sufficient to explain the standard test data available 
at that time (Dicke 1957). 

The Yilmaz theory is a candidate for the generalized - anisotropic - version of PV. 
As remarked in the introduction, the Yilmaz theory shares with PV the electromagnetic 
motivation, and, in the case that only Too is non-zero, the same form for the metric. This 
suggests that PV may turn out to be a truncated version of the Yilmaz theory, in the 
same sense that a theory of electromagnetic interaction would be truncated by keeping 
but ignoring A. Since the Yilmaz theory has the full compliment of 10 functional degrees of 
freedom in the metric, it has the potential to resolve the problem of the reduced radiation rate 
suffered by PV. To date there have been no published studies of the cosmological predictions 
of the Yilmaz theory, so it would be interesting to pursue this matter further in a future 
report. Since the Friedmann equation is predicated on isotropy and homogeneity on the 
cosmological scale, it is conceivable that such a development will affect only the compatibility 
with the binary pulsar data through the polarization states available to radiation, without 
affecting the cosmological predictions of PV reported above. However, no conclusions are 
possible based on these very broad arguments. 



5.2. Vacuum contribution 

It is recalled that for the Einstein equation with a Cosmological term to retain its tensor 
form, it is necessary that the vacuum action must have the form constant x y/—g. With this, 
the FLRW flat-space metric (i.e. with line element equation (13)) gives that the vacuum 
contribution appears in the Lagrangian density as constant x a (t) 2 (equation (6)). The PV 
field K, however, is a scalar, and therefore its (Euler-Lagrange) equation of motion is a scalar. 
It follows that there is no such constraint on the vacuum contribution to the PV action. In 
principle therefore, other vacua may be contemplated, with potentially radical consequences 
for the resulting cosmology. In the absence of a solid mathematical basis for a particular 
form for the vacuum term, some other argument in favor of a particular alternative would 
have to be provided in order to proceed further. It is noteworthy though that PV cosmology 
is in closest agreement with observation when the vacuum term is close to zero. Given this, 
and the arbitrariness of the function dependence of the vacuum term on the scale factor, it 
could be argued that the vacuum term be expunged from PV altogether. 
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Alternatively, one might argue that with no constraints on the vacuum dependence on 
the scale factor there is some justification for regarding the total energy term in equation 
(16)) as a novel, negative-energy vacuum term. That is, equation (6) could be replaced with 

L vac -> - J d 4 xa 2 p v (a); p v (a) = -p v /a 2 (71) 
where now p v is a positive constant, which would then generate the Friedmann equation 

Hq \dr J a a 2 

The advantage of this view is that there is now no arbitrary overall energy - the total energy 
is now zero. Of course, the arbitrariness in f2 s has just been subsumed into the arbitrariness 
of the vacuum. 



5.3. PV with a spatially- varying 7^-field 

Recalling that the PV action is as if the matter and fields are in a curved space-time with 
line element (3), (of which the flat-space FLRW line-element of equation (8) is obviously, 
therefore, a special case), it is observed that Dicke's scalar PV is committed to an isotropic 
space. In the case that K is just a function of time - as it is here - then equation (8) is the 
FLRW flat space metric (with proper time r then given by dr — dt j \J~K (£)= dt/a (t) ). It 
is obvious from the form equation (8) that the other FLRW space-geometries (spherical and 
hyperbolic) are not accommodated PV. The line elements corresponding to those geometries 
were first derived independent of GR (see for example Tolman (1987)) based upon Milne's 
Cosmological Principle, which, if regarded as an axiom, dictates that PV cosmology entertain 
only the flat-space metric. 

It is noteworthy that the specific form of the three FLRW line elements (plus, of course, 
the Friedmann equation) can be generated from 'within' GR, i.e. without direct appeal to 
the Cosmological Principle, and without enforcing any symmetry by hand. For example, it 
sufficient (in order to generate all FLRW three metrics), that that the matter is required to 
be static and uniformly distributed (in whatever space it exists), and that the metric can be 
written in the form 

ds 2 = c 2 dt 2 - f (t, x) dx 2 (73) 

Given that this derivation of the FLRW metrics is, in some sense, achieved more 'directly' 
by the theory of gravity (rather than by a symmetry principle), it could be argued that this 
is the physically superior of the two routes to the metric. From that standpoint, a similar 
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procedure should be used to decide on the cosmological metrics appropriate to PV. To that 
end, alternatives to the flat space assumption might arise from an 'effective' cosmologically- 
appropriate version of the PV action, equation (4), that does not, a priori, exclude spatial 
variation of K: 



= -/"**{ 
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Variation of K gives the Euler equation 

72^ /V7^\ 2 Kd 2 K 
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(74) 



dt 2 



Pv 



Pm Pr 



(75) 



where the pi are constant. In principle, given plausible boundary conditions for, say, K 
and dK/dt everywhere at t — 0, solutions to equation (75) should provide cosmologically 
acceptable functions of K that vary in space and time and which will generate PV analogues 
of the FLRW metrics. 



5.4. Relation between spatially- varying 7^-field and vacuum and overall 

energy densities 

In the Hamiltonian corresponding to the action in equation (74) the spatial gradient 
term contributes positively, with the same sign as the matter and radiation terms (the sign 
of the vacuum term is not fixed). Consequently any spatial variation in the metric increases 
the overall system energy. This is to be contrasted with the effect of the spherical and 
hyperbolic FLRW metrics in GR, which contribute terms of opposite signs to the Friedmann 
equation for the scale factor a(t). Since agreement with observation in a flat-space zero 
vacuum cosmology already requires that PV have non-zero positive energy (i.e. ~ 0.3 
identified above), it follows that the introduction of space-dependent forms for K cannot be 
used to offset this requirement. In other words, PV without a contribution from the vacuum 
cannot have zero overall energy in any space-time. 



5.5. Possibility of a radiation era 

Analysis above showed that flat-space PV has no baryon, lepton, or radiation eras. 
Though the baryon and lepton eras were shown to be orders of magnitude beyond the reach 
of PV, the same was not true of the radiation era, which was discounted on the basis of an 
inferred age of 6.6 Gyr or less. It is conceivable that a spatially varying if-field satisfying 
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equation (75) with appropriately-chosen boundary conditions could satisfy the observational 
constraints on age etc, and bring a radiation era within the scope of PV by a small adjustment 
of the metric at very early times (yet without significantly changing the later evolution). 

5.6. Gravitational redshift and Yilmaz stars 

With reference to equation (75), the vacuum behaviour of the K -field near a massive 
object (ignoring Cosmological influences) is given by 

V 2 K/K -(VK/Kf = 0. (76) 

Setting ip = log if gives the Laplace equation, V 2 ?/> = 0, from which it is inferred that 

oo / 
1=0 m=-l 

where the Yi m (9, (f>) are spherical harmonics, and where ai m and b\ m are arbitrary constants 
to be determined by the boundary conditions. But far from the source of the field K must 
return to a constant value, so all the ai m and b 0m must be zero, and therefore 

oo I 

k = n n ex p ( bimYim ( e > ^/ ri ) ( ?s ) 

1=1 m=—l 

where the remaining b\ m can be determined by the conditions at the surface of the source. 
The important feature of equation (78) is that, regardless of the mass distribution, K is 
everywhere finite outside the mass. Recalling the PV line element, (3), it is seen that the 
coefficient of the dt 2 term in the line element is never zero, from which it is deduced that PV 
does not admit event horizons outside the matter distribution; PV predicts that light from 
the surface of any star, however massive, cannot be infinitely redshifted, i.e., a star can grow 
to an arbitrary density without disappearing from view. Hence it is concluded that PV does 
not admit black holes. This statement does not address the issue of stability; properly applied 
with an appropriate equation of state, PV may or may not concur with GR predictions for 
instability and collapse. At this point the only certain (but important) difference is that any 
collapsed star will nonetheless remain visible, albeit with an arbitrarily large redshift. As a 
result of this difference with GR, a more thorough appraisal of the validity of PV might need 
to take into account the possibility that GR-designated cosmological redshifts may have a 
significant gravitational component. Since, as pointed out in section 5.1, the exponential 
form of the vacuum metric, equation (78), is common to PV and the Yilmaz theory, the lack 
of black holes is a prediction common to both theories. 
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5.7. Relation to the quasi-steady-state cosmology 

Several authors have presented arguments in favor of a quasi-steady-state cosmology 
(QSSC) that has no big bang (Burbidge 2001; Narlikar, Vishwakarma & Burbidge 2002; 
Volkov, Izmest'ev & Skrotskii 1971). A presentation of both sides of the QSSC-GR debate 
may be found in Burbidge, Hoyle & Narlikar (1999) and Albrecht (1999), respectively. Lack- 
ing BBN, they have argued that observed abundances of the light elements are the result of 
stellar rather than big bang nucleosynthesis - a feature of potential relevance to PV. Also in 
common with PV, QSSC also lacks an era of radiation domination. QSSC claims that the 
CMB is thermalized starlight, thermalisation occurring principally in the most contracted 
phase of an oscillating cosmology. This component of the QSSC hypothesis is potentially 
relevant exclusively to the oscillating PV cosmologies, including the particular example dis- 
cussed above and shown in Fig. 4 (which happens to have the same oscillation period, of 
around 100 Gyr, as that preferred by QSSC). The degree of thermalisation of light from a 
previous cycle will depend on the rate and depth of the contraction, the details of which 
differ in QSSC and PV. Further investigation would be necessary to test the viability of this 
idea within PV. 
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